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Abstract: With an appropriate choice of parameters, a higher derivative theory of gravity 
can describe a normal massive sector and a ghost massless sector. We show that, when defined 
on an asymptotically de Sitter spacetime with Dirichlet boundary conditions, such a higher 
derivative gravity can provide a framework for a unitary theory of massive gravity in four 
spacetime dimensions. The resulting theory is free not only of higher derivative ghosts but 
also of the Boulware-Deser mode. 



Contents 



1 Introduction 



1 



2 Higher derivative ghosts 

2.1 Scalar toy model 

2.2 Massive gravity with boundary conditions: Linearized analysis 



2 



3 
4 



3 Boulware-Deser ghost: Full ADM analysis 



7 



4 Discussion 



9 



1 Introduction 

The seemingly simple task of giving mass to a graviton or a spin-2 field has proved to be 
rather challenging. To obtain a linearized theory of massive gravity, we add a mass term 
to the linearized Einstein-Hilbert action. Out of the ten degrees of freedom (DOFs) of the 
4D metric perturbations, four are removed by requiring covariant conservation of energy and 
momentum. But one of the remaining six DOFs being a ghost, we are forced to choose a 
special form of the mass term, to get the so-called Pauli-Fierz (PF) theory [1] where an extra 
on-shell symmetry kills the sixth mode. Then, 

• A massive graviton has more DOFs than the massless one of general relativity (GR), and 
at the linearized level its coupling to matter is finitely different from that of GR even 
for arbitrarily small graviton mass. This discrepancy, dubbed as van Dam-Veltman- 
Zakharov (vDVZ) discontinuity, leads to conclude that massive gravity would be exper- 
imentally disfavored [2]. 

• One of the five DOFs of the PF theory can get strongly coupled at too low energy 
scales [3]. While this affects the predictivity of the theory, it also signifies that the 
vDVZ discontinuity - that was derived assuming the linearized theory is valid at short 
scales - should not worry us too much. 

• Upon turning on generic higher order interactions, the extra symmetry of the PF theory 
gets broken, and the sixth mode or the Boulware-Deser (BD) ghost reappears [4]. 

• The structure or tuning of the Pauli-Fierz mass term is not radiatively stable, and loop 
corrections from generic quantum interactions will destroy it [5]. 

• It is highly non-trivial to construct a consistent non-linear completion of the Pauli-Fierz 
theory. 
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The study of massive gravity is not confined to the somewhat academic chahenges hsted 
above: interest in this field has been resurrected by the advent of theories of modified gravity 
as an alternative way of tackling the cosmological constant problem - in particular, some of the 
popular higher dimensional modified gravity models are, from the effective 4D theory point 
of view, PF theories. Therefore the construction of a theoretically and phenomenologically 
consistent model of massive gravity could affect the discussion about the accelerated expansion 
of the Universe. 

Recently there have been two major developments in the quest for a consistent theory of 
massive gravity. First, ref. [6] formulated a non- linear completion of massive gravity without 
the BD ghost by introducing an auxiliary metric and an elaborately designed interaction 
potential between the physical metric and the auxiliary one. The second [7] employed higher 
derivative gravity (HDG) as a framework of a consistent massive gravity in 3D: HDG once 
waxed due to its capability of addressing renormalizable quantum gravity [8] , but then waned 
because it generically contains extra ghost DOFs that destroy the unitarity of the theory. 
But there being no propagating massless gravitational DOFs in 3D, HDG in 3D has only the 
massive sector, which can be chosen to be non-ghost-like. 

Extending the scheme of [7] to 4D is not trivial, because in 4D the massless sector still 
has two DOFs; if we choose one of the sectors to be normal, the other is bound to contain 
ghosts. Then there came [9], whose author took the conformal gravity of (Weyl tensor)^ 
which, due to its higher derivative nature, has both a massless and a massive tensor with 
opposite sign kinetic terms, and put it in a maximally symmetric spacetime. By imposing 
appropriate boundary conditions to truncate massive ghost modes, [9] showed that conformal 
gravity can reproduce GR at the classical level. One can expect that combining the ideas 
of [7] and [9] might provide a scheme to obtain a unitary theory of massive gravity in 4D.^ 
In the next section we will show that this scheme does indeed work for HDG on a de Sitter 
(dS) background. Then in §3, we check that HDG can be free of one of the most notorious 
problems of massive gravity, the BD mode. We briefly discuss more of the possible issues in 
§4. 

2 Higher derivative ghosts 

In order to make sense of a theory defined with higher time derivatives, we need a way to 
project out the unwanted ghost DOFs. We will show that for a theory defined on a dS 
background, imposing Dirichlet boundary conditions on the dS boundary provides such a 
projection mechanism. 

^This path is explored in [10] to construct a theory of a partially massless spin-2 field in a conformal gravity 
setting. 
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2.1 Scalar toy model 

Let us first warm up by considering a higher derivative action for a scalar field ^, which we 
write in the form 

S^ = \j d^x^{in-i>f + m^d^^V>'^Y (2.1) 
in a dS background, ds^ = ^_Ht)'^ (— dr^ + dx^). By introducing 



m 

we can rewrite it as 



— , ^P = ^ + ^, (2.2) 



5^ = ^ / d^xV^{(VV^)2 - (V<A)2 -m'<p']. (2.3) 

This expression shows that the action (2.1) does actually describe two scalar degrees of 

freedom, and that, for m? > 0, the massless mode has a "wrong sign" kinetic term, i.e., is a 
ghost field. Now we solve the ip- and (/>-equations of motion (EOMs): 

□ ^ = 0, (□-m^)(/) = 0, (2.4) 

or 

^"-l^' + kH = o, 4>"-l4>' + {k' + -^)4> = o, (2.5) 

where " implies a Fourier transform along the three spatial dimensions and ' = The 
solutions are 

= (-r)3/2{ai Hff^{-kT) + 6i Hfli-kr)] , (2.6) 

^ = (-r)3/2{a2i7()) (-^^) + b2Hf^^{-kr)] , (2.7) 



with 1/ = y^9 — Am? /H'^. We will use as mode functions the positive frequency modes at 
small scales (|/cr| » 1), so that we set 6i = 62 = 0. 

To obtain a unitary theory we need to project out the massless mode tp. To this end we 
observe that, since the behavior of the solutions as r — t- 0~ is 

^^_i(_,)(3-.)/2(|y"=r{^)a.. (2.9) 

imposing the Dirichlet boundary conditions, "0^=0 = and </>|r=o = 0, forces ai = 0, while 
the (^-mode survives as long as Re(z/) < 3, i.e., m? > 0. In terms of the original field ^, 
this amounts to imposing the Dirichlet boundary condition lim.^_^Q- (^_ij^yi ^ = 0. There- 
fore, although we started with a higher derivative theory, by imposing appropriate boundary 
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conditions and by choosing the overall sign of the action (2.1) such that the massive field (f) is 
normal, we can obtain a unitary theory of a single massive field. Once we have selected the 
linearized modes of the theory this way, the effect of sources or of self-interactions of the ^ 
field can be accounted for in a perturbative sense, and the theory is guaranteed to be ghost- 
free at tree level. Of course, one has to make sure that the sources for ^ allow consistent 
imposition of the boundary conditions. 

2.2 Massive gravity with boundary conditions: Linearized analysis 

The most general action for an "i?^"-HDG in 4D contains a cosmological constant, i?, B? and 
R^,yR^"' {R^.^pxR'^'P^ can always be traded for R'^ and R^.uR^", 

since in 4D the Gauss-Bonnet 

invariant is a total derivative) and can be shown [11] to describe a massless tensor, a massive 
tensor and a massive scalar. We will assume a special relation between the coefficients of 
Rf^uR^'^ and R?, with which the scalar mode gets infinitely heavy and is removed from the 
theory.^ We will discuss in §4 our expectations in the case where R^yR^'^ and R? carry 
arbitrary coefficients. 

Based on the considerations above, we choose the action of our system to be 

I- |- , 1 \ 1 

S = ^ d'^Xy/^^aR- A + a [R^uR^'' - ^-R^jJ + boundary terms , (2.10) 

where a = +1 or —1, which will be determined by requiring the massive sector to be normal. 
We perform a perturbative expansion of (2.10) on top of the de Sitter background 

1 



dg' = (_ j^^)2 (^M^ + V) dx^dx" . (2.11) 
At 0{h), tadpole cancellation fixes A to be 

A = 6aH'^. (2.12) 

In order to count the number of the physical DOFs, we decompose /i^j^ into different 
helicity modes: 

/loo = -2 </> , hoi = + diB , 

hij = hj^ + d^(f + djij -25,j^ + 2 d^djE , (2.13) 

where dih^^ = 0, ^ijhj^ = 0, dihj^ = and di^J = 0. By construction, there is no mixing 
between different helicities at the quadratic level, and we can consider the helicity~2, -1 and 
-0 sectors separately. 

The helicity-2 sector is simple: 



2 



(2.14) 



^In 3D, the same consideration would fix tlie ratio of tfie coefficients of R^^R'^^ and R? to be —3/8, which 
is the choice leading to the New Massive Gravity of [7]. 
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with A = df. Introducing the auxihary field fij to ehminate higher r-derivatives turns (2.14) 
into 



(2) 

hol-2 



Ml 



d^x 



a-2H^a 



+ —(f'.. /iTT/ , f. . A /,TT^ _ fZ 



(2.15) 



and diagonahzing it gives 

M2 ^2 ^a-2H^a) 



(2) 

hcl-2 



d^X 



r 1 



1 



H \2 



2 - wp-jH"' 



J- 



(2.16) 



where hT^ = -Ht-i^ {X^j + Yij), fi, = -Ht- ^(2- Yij, = 2 - ^ a.nd is 

an arbitrary normahzation constant. To determine which modes are ghosts, we employ the 



Hamiltonian formulation with S*, 



(2) 

hcl-2 ~ 2 2 



(2) 

hel-2 



Ml ^2 



2 2 
+PyY' 



J d^x[PxX'-{ 



p2 



2{a-2H^a) 



hcl-2- 



2(a-2H^a) 



2 Itt2 



+ ^ y A + ^ 



}], (2.17) 



5C 



(2) 



where Px.y = 5[x,y) ^^^'^ "^3 indices are suppressed. Therefore in order for the Hamiltonian 
of the massive helicity-2 modes to be positive definite, we have to require 



a -2H^ a <0. 



(2.18) 



Notice that this implies the massless helicity-2 modes are ghost-like. 

The helicity~l sector, once the gauge invariant variable is identified, is even simpler. 
With Vi = ^/^{hji - (['), we have 



5, 



(2) 



Ml 



a d X 



-(.D^ + -.,.(a+ ^/ 



(2.19) 



-"hel-l ~ 2 

From the Hamiltonian analysis similar to that of the helicity-2 modes, a must satisfy 

Q>0. (2.20) 
The helicity-0 sector requires some work. Employing the gauge invariant variables 

$ = (/,-tV', B = A{B - E' + Tip), (2.21) 
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the helicity-0 part of (2.10) is^ 



'hel-O 



(2) 

Note that $ is non-dynamical. With 5'j^pj_o 
formulation gives 



*4 f^4^r(2) 



5, 



(2) 

hel-O 



(2.22) 

2 / d^x£^g{_Q, writing in Hamiltonian 
2a 



> a 



if2 r4 



$(2ri3-3$) 



Since $ is a constraint, we integrate it out to get 

3 A2 



(2) 



d^x 



-^hel-O 2 

Upon a canonical transformation, 

2a 



24 o- 



Pb 



V3H- 



2a 



we finally obtain 



^(2) 

'hel-O 



Mi 



p2 



-a) 



which dictates a to be negative, i.e., 



2{-a) 



a 



2 ^ 2^ 
A + ^ 



Q. 



q} 



-1. 



(2.23) 



(2.24) 



(2.25) 



(2.26) 



(2.27) 



Let us then combine the results from all three sectors. The positivity conditions for 
three sectors, (2.18), (2.20) and (2.27), can be simultaneously satisfied, and we choose // = 
2 _g'^^Q+ T ' -^2 = Y, Fi = ^/av and Fq = Q, to write the full quadratic action as 

. f x; . i.„(A . ^).,.} - 1.^ - ix(A . 



with 



2H^ + 



1 



a 



(2.28) 



(2.29) 



That is, there are a normal massive spin-2 field with 5 DOFs and a ghost massless helicity-2 
mode with 2 DOFs. Note that the Higuchi ghost [12] is not present because > 2H^. 



^As mentioned earlier, instead of (2.10) if we had started with 

J d^x^i^aR - A + aR^^R"-" - I3R'') 



with a 7^ 3/3, there would have been two propagating scalar DOFs. In the current formulation, this would 
have been manifested by the presence of a higher r-derivative term: (a — 3/3) (i/)")^. 



- 6 - 



The solutions of the equations of motion for the Fourier modes Fn and X are 

Fn = ^{anHl^/^i-kr) + bnHlJ^{-kr)} , X = ^ {c H^^/^{-kT) + dH^f^i-kT)} , 

(2.30) 

with = Y^9 — Am? /H"^ = -y/l — 4/(a H"^). Requiring the standard short scale behavior 
immediately sets 6„ = d = 0. 

As we did for the scalar toy model, we now use boundary conditions to eliminate the 
unwanted ghost. ^ Since we impose boundary conditions on the metric perturbation h^i, and 
the leading r-dependences for small |t| are F2, X ~ '^~^^Jj^^ ~ h^^ and Fq ~ B, we need 
to know the behavior of TF2, tX, Fi and as r — )■ 0~: 

A ^ -l(-r)<-")/^(f)"''r(^)«„ (2^32) 
^o^4(-.)"-'"'^(|)''"r(^)ao, (2.33) 



Therefore we see that, if we impose the Dirichlet boundary condition on h^i, at r = 0~, the 
five massive modes can survive because Re(i/) < 1, but the massless ones are projected out, 
and we are left with a unitary massive spin-2 system. 



3 Boulware-Deser ghost: Full ADM analysis 

In addition to the higher derivative ghost, HDG, like any theory of massive gravity, might be 
haunted by the BD ghost. Since the BD mode shows up at the nonlinear level, the best way 
to see if it is there is to employ the full Arnowitt-Deser-Misner (ADM) formulation [14] and 
count the number of Hamiltonian constraints. 

With the introduction of an auxiliary metric, f^,^, (2.10) can be written as 

S = ^l d'xV^[aR + U4R^'^-\g^''R)-^if,.F''-f)}. (3.1) 

For ADM treatment, we first decompose the metrics into ADM variables: 

_ (-N'^ + NiN' N,\ _ (-2Nn + 2Nin' + ^ijN^N^ m + ^ikN^" 

*See, e.g., [13] for issues regarding imposing future boundary conditions in dS spacetime and their possible 
resolutions. 
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to obtain 



Ml 



4a V^^^'^ ^ N 



+2Kk-2 Kij k'^ + ^ Vin* - K^^ Vin^) } , (3.3) 



where 



Kij = ^(-r^i + ViNj + VjNi) , (3.4) 

kij = ^{-iij + N'' Vk-f^j + 7fci iV^ + jkj ViiV'^) . (3.5) 

Three-dimensional indices are raised and lowered by Tij , so that K = F*-' Kij , k = T^^ k^j and 
7 = r*-' -jij. By defining the canonical momenta conjugate to Tij and 'jij by 



we can put (3.3) into the canonical form 

s = ^ y d4x(p^^' r,, + j,j - n) , (3.7) 

with 

n = -2ni VjP'^ + {P'^ Vaij - 2 Vj{P'^J^k) - 2 ViUl] 

32 



^ { - ^ - i) ^^^^ + - + - 4^ {^^^ 

{2 n,, - p n + (a + ^ - I) (p^^- p,,. - + 2 p,, p^'^ 7I - 7., } • (3.8) 



After integrating out the constraint rii, we finally obtain 



n = Nno + NiUi + nno, (3.9) 

where 

= -^{Pn - 2P^^ni, + (a - I) - P^^P,,) + PP^^7,j - 2Pi,P^Sl^ + 2aViP^*^VjP^ } 

_ Vf { [a - 1) (3)P + A7 - V,V,7'^' + (')ii^j7*^' - ^(7*^^' - 7')} , (3.10) 
Hi = P'^'Va.k - 2V,(P^Sifc) - 2V,n^' , (3.11) 
7io = Vf((3)p-l7)-^(P-P.-^). (3.12) 
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The lapse, N, and the shift, Ni, for g^^, are Lagrange multipliers, as they should be: the 
resulting constraints immediately remove four DOFs, and additional four will be eliminated 
with secondary constraints. What is important to note is that the lapse, n, for ffj_u is also a 
Lagrange multiplier, so that two more DOFs will be killed by the corresponding primary and 
secondary constraints. Therefore, DOF counting goes as 

6[Tij] + 6[Pi^] + 6[7,j] + Q[U,j] - 2 X (1 + 3 + l)[iV, iV„ n] 
= (5 + 5) [massive graviton] + (2 + 2) [massless graviton] , 

showing that, even at the full nonlinear level, the action (2.10) describes 5+2 propagating 
degrees of freedom and there is no 6th massive mode, i.e., the BD ghost. Then imposing 
boundary conditions will remove the massless modes in the end. 

4 Discussion 

In this paper we have described a way of realizing a non-linear completion of PF massive 
gravity without ghosts. Yet, there are more items in the to-be-resolved list presented in the 
introduction. Let us briefly discuss a couple of them. 

The vDVZ discontinuity is not supposed to be an issue in the present set-up, because the 
discontinuity disappears on (A)dS backgrounds as m?/H'^ — )• [15]. However, in the case of 
dS background we have considered here the range < m? < 2 is forbidden, since massive 
graviton would be a ghost in this regime; the conditions (2.27) a = —I and (2.20) a > were 
imposed precisely to avoid the appearance of such a ghost. Therefore, the limit m — t- does 
not belong to the allowed parameter space for our scenario. 

It is worth noting that HDG ("curvature squared gravity", to be precise) shares a set 
of solutions with GR. In particular, in the presence of a spherically symmetric source, the 
Schwarzschild-dS (SdS) metric solves the EOMs of HDG and satisfies the Dirichlet boundary 
condition on dS boundary, which seems enough to make SdS a viable solution for our system. 
But it turns out not to be the case. In fact, our EOMs can be schematically written as 

-Gi,^ + -gf,^ + a£^^ = j^, (4.1) 

where Gfj^u is the Einstein tensor, the contribution from the curvature-squared terms and 
Tfj^u the stress energy tensor for a spherically symmetric source. For solutions of pure GR, 
£^1, identically vanishes, and we are left with an Einstein-like equation with an overall minus 
sign on the left hand side, which results in a negative Schwarzschild radius. If we do not allow 
for naked singularities in the theory, SdS is not a valid solution of our system. It would be 
interesting to find a physical, asymptotically dS, spherically symmetric solution to (4.1) and 
to analyze its behavior at various scales. 

As for the stability/robustness of the model, we stress that we are considering it only at 
the classical level (this is not more restrictive than what happens in the other massive gravity 
theories discussed in the introduction). In a general quantum mechanical theory, indeed, we 



- 9 - 



would expect R^, R^, ■ ■ ■ corrections to be generated and to affect with 0(1) corrections the 
properties of the massive mode. A related question concerns the robustness of the theory with 
respect to the tuning of the ratio, — |, of the coefficients of R^^R^'^ and R^ in the lagrangian, 
i.e., whether we can still have a unitary theory even if this tuning is broken. When we start 
with J d'^Xy/—g(^a R — A + a R^u R^^ — /3 i?^) , the theory contains an "eight mode" with 
mass proportional to (a — 3/3)~^ [11] that, for our choice a = —1, should behave as a ghost. 
Now, in analogy with the analysis of section 2.2, one may speculate that the scalar ghost 
could be truncated by the Dirichlet boundary conditions - at least in part of the parameter 
space spanned by a and /3. It might be worthwhile to confirm this speculation with a full, 
explicit investigation. 
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